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Using the stabilized jellium model with self-compression, we have calculated the dissociation 
energies and the barrier heights for the binary fragmentation of charged silver clusters. At each step 
of calculations, we have used the relaxed-state sizes and energies of the clusters. The results for the 
doubly charged Ag clusters predict a critical size, at which evaporation dominates the fission, in good 
agreement with the experiment. Comparing the dissociation energies and the fission barrier heights 
with the experimental ones, we conclude that in the experiments the fragmentation occurs before the 
full structural relaxation expected after the ionization of the cluster. In the decays of Ag^ 1 " clusters, 
the results predict that the charge-symmetric fission processes are dominant for smaller clusters, 
and the charge-asymmetric fission processes become dominant for sufficiently larger clusters. 
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I. INTRODUCTION 



The fragmentation of ionic metal chistersi'i as well as other properties of metallic clusters have been extensively 
studied using the jellium model (JM).H~0 In this model, the discrete ions are replaced by a uniform positive charge 
background of density n — S/Anr^ in which r s is the bulk value of the Wigner-Seitz (WS) radius of the valence 
electrons of the metal. The simplest geometry for the positive background is spherical which is appropriate for closed- 
shell clusters or large enough clusters in which the Jahn-TellerQ deformation has negligible contribution. However, 
using this simple spherical JM, a iot of information on the properties of metal clusters has been obtained. A refined 
version of the JM, the stabjJizedQQ jellium model (SJM), which was introduced by Perdew et. al. in 1990, has 
improved some drawbacksBEJ of the JM (For a-pecent review on SJM,-see Ref . Jl3| ) . Inr-Mceiit years, the SJM has-been 
used to predict the properties of bulk metalsQO, metal surfaces! 1 To, metal clustersEjtll and metallic voids.EJ The 
fragmentation of charged metal clusters has been also studiedE3E3 by Vieira et. al. using the SJM. However, since 
the surface effects have a large contribution in the energetics and sizes of small clusters, and also since in a fission 
process the competition between the surface tension and coulomb repulsion leads to the existence of a barrier, a more 
sophesticated useOO of the SJM is needed to predict the correct energetics of the clusters and the barrier heights 
(BH) in the study of the fragmentation processes. This method, which is, called SJM with self-compression (SJM-SC), 
has been used to predict the equilibrium sizes and energies-of charged^ or spin-polarizedEZI metal clusters as well as 
the calculation of chemical potentials of metallic clusters.EHl The SJM-SC has been also used by Sarria et. alJtn to 
calculate the surface energies and the work functions of metals. In contrast to the JM and the SJM in which the r s 
value is borrowed from the bulk system, in the SJM-SC, the density parameter r s of the jellium sphere is assumed to 
be a free parameter which can be adjusted in such a way that a cluster with a given number of electrons an d spe cific 
electronic configuration achieves its equilibrium state. The SJM-SC calculations on neutral metal clusteralaO has 
shown that the equilibrium r s value of the jellium sphere is less than the bulk value and tends to its bulk value for 
infinitely large cluster. This phenomenon is called self-compression which is due to the dominant effect of surface 
tension in small metal clusters. However, it has been shown thatEil charging a small metal cluster can result in an 
equilibrium r s value which is larger than the bulk value-pThis effect is called self-expansion. The self-expansion has 
been also predicted for highly polarized metal clustersE3E3. These two effects have different origins. In the former, 
the repulsive coulomb force dominates the surface tension whereas, in the latter, the Pauli force is responsible for the 
self-expansion. 

In this work, using SJM-SC, we have studied the binary decay processes of positively charged Ag clusters containing 
up to 100 atoms in all possible channels. We have considered the following possible decay processes for singly ionized 
Ag clusters 

Ag^->Ag&_ p + Ag°, p = l,2,-..,JV-2. (1) 

For doubly charged clusters, the decays can proceed via two different processes. The first one is the evaporation 
process 

Ag^-AgJ+p + Ag", p=l,2,---,N-3 (2) 

and the second one is fission into two charged products 

Ag^+^AgJ+p + Ag^, p = 2,3,.-.,[JV/2] (3) 

In general, for the binary decay of Z-ply charged (Z is a positive integer) cluster, we have 

Ag^ - Ag^ z ; +Ag^, « 1 = 0,l,-..,[Z/2]; p = z x + 1, • • • ,N - Z + z x - 1. (4) 

For an even value of Z with z\ = Z/2, the range of p reduces to p = z\ + 1, ■ • • , [N/2]. The processes for which z\ = 
( i.e., one of the fragments is neutral), are called evaporation processes and others ( both fragments are charged) 
are fission processes. In evaporation processes, the negativity of the difference between total energies before and 
after fragmentation is sufficient to have a spontaneous decay. However, in fission processes a negative value for the 
difference energy is not sufficient for the fission of the parent cluster. This is because, the competition between the 
short-range surface tension and the long-range repulsive coulomb force may give rise to a fission barrier ( i.e., one 
should supply energy to overcome the barrier). 

The organization of this paper is as follows. In section || we explain the method of calculating the total energies and, 
fission barriers. To obtain the total energy of a given cluster, we solve the self-consistent Kohn-Sham (KS) equationsc.il 
the density functional theoryo (DFT) with local spin density approximation (LSDA) for the exchange-correlatio 



m 



(XC) functional. To calculate the fission barrier, we use the two-touching-spheres model for the saddle configuration 



In section [II, we discuss the results, and finally, we conclude this work in section [V 
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II. CALCULATIONAL SCHEME 



A. Total energy of a cluster 



In the context of the j=UM, the average energy per valence electron in the bulk with density parameter r s and 
polarization ( is given byE!3 



e(r s , C, r c) = t s {r s ,C) + £xc(r s ,C) + w R (r s ,r c ) + e M (r s ), 



where 



t s (r s ,Q = ^ (1 + C) 5/3 + (1 - 



\5/3 



(5) 

(6) 



£xc (r s ,0 = + C) 4/3 + (1 - C) 4/3 1 + e c (r s , () 



(7) 
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All equations throughout this paper are expressed in Rydberg atomic units. Here t s and e xc are the mean noninter- 
acting kinetic enecgy and the exchange-correlation energy per particle, respectively. For e c we use the Perdew-Wang 
parametrization.E3 For a z-valent metal the average Madelung energy, em, is defined as £m = — 9z/5ro, in which ro is 
the radius of the WS sphere, r — z 1 / ? 'r s . In Eq.([3|), £ — (n^ — n;)/(ri| + in which and are the spin densities 
of the homogeneous system with total density n_^j= nf + ri[. The quantity wr is the average value (over the WS cell) 
of the repulsive part of the Ashcroft empty coreo pseudopotential, 
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(9) 



and is given by wr = 3r^/r^ where, z is the valence of the atom, 9{x) is the ordinary step function which assumes 
the value of unity for positive arguments, and zero for negative values. 

The core radius is fixed to the bulk value, rf , by setting the pressure of the unpolarized bulk system equal to zero 
at the observed equilibrium density n = 3/47r[ff (0)] 3 : 



_d_ 



e(r s ,0,r c ) 



0. 



(10) 



r s —rf(0),r c —r^ 1 



Here, ff(0) = rf(C = 0) is the observed equilibrium density parameter for the unpolarized bulk system, and takes 
the value of 3.02 for Ag. The derivative is taken at fixed r c , and the solution of the above equation gives as a 
function of ff (0) 



'(0)] 
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d 
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(11) 



; (o). 



The SJM energy for a spin-polarized system with boundary surface is given by! 

ikiM [nt,n;,n+] = £,tm [n-i, n;, n+] + (e M (r s ) + w R (r s , rf )) J drn+(r) 
+<«»)ws(r.,rf) J dre(r) [n(r) - n+(r)] , 

where 

Ej M [nt , n;, n + ] = T s [n T , n{\ + E xc [ra T , n;] 

+ - / dr (j>([n,n+];r) [n(r)-ra+(r)] 



(12) 



(13) 
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^([n,n + ];r)=2/dr'M^±^. (14) 



and 



Here, n = n-\ + and n + is the jellium density. 0(r) takes the value of unity inside the jellium background and 
zero, outside. The first and second terms in the right hand side of Eq.(^) are the non-interacting kinetic energy and 
the exchange-correlation energy, and the last term is the Coulomb interaction energy of the system. The quantity 
(<5i>)ws is the average of the difference potential over the Wigner-Seitz cell and the difference potential, Sv, is defined 
as the difference between the pscudopotcntial of a lattice of ions and the electrostatic potential of the jellium positive 
background. The effective potential, used in the self-consistent KS equations, is obtained by taking the variational 
derivative of the SJM energy functional with respect to the spin densities as 

<SS ([»T.»i. «+] 5 r ) = s^)^ Esm ~ Ts) 

= <P ([n, n+] ; r) + < c ([n T ,nJ ; r) + Q(r)(5v) ws (r s , rf), (15) 
where a = \, \ . By solving the KS equations 

(V 2 + < // (r))^(r) = £ f^(r), * =U, (16) 

n(v) = M^), (17) 
r* ff (r) = ICWI 2 , (18) 

i(occ) 

and finding the self-consistent values for ef and </>f , one obtains the total energy. 
In our spherical JM, we have 

n +(r) = ^0(R-r) (19) 

in which R = (zTV) 1 / 3 ?^ is the radius of the jellium sphere, and n(r) denotes the electron density at point r in space. 
Using the Eq. (21) of Ref. §, this average value is given by 

Wws(r,,r?) = -Ui--— . (20 ) 

Applying Eq. ( |l^ ) to a metal cluster which contains spin-up, N± spin-down and N (= + Ni) total electrons 
in the ground state, the SJM energy becomes a function of N, r s , and rf . The equilibrium density parameter, f s (N), 
for a cluster in the ground state electronic configuration, is the solution of the equation 



= 0. (21) 

r s =f s (N) 



Here, the derivative is taken at fixed values ofjjy and rf . For an iV-electron cluster in its ground state electronic 
configuration, we have solved the KS equationscj self-consistently for various r s values and obtained the equilibrium 
density parameter, f s (N), and its corresponding energy, E(N) = E(N,f s (N),rf). 

B. Dissociation energy and fission barrier 

The dissociation energy (DE) for the general binary decay process (||), defined as the difference in the sum of total 
energies of the products and the total energy of the parent cluster, is given by 

D Z Z1 (N, p) = {E z N -_!-; +E^)-E Z , (22) 
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In evaporation processes (zi = 0), a negative value for the DE implies that the parent cluster is unstable against 
that particular decay channel and therefore, the fragmentation is spontaneous. On the other hand, a positive DE in 
a particular decay channel means that the parent cluster is stable against the decay in that particular channel. That 
is, one should somehow supply energy to the system to induce the fragmentation. For processes (|l|) and (||) the DE 
becomes 



and 



Dl + (N,p) = (Ej+ + El) E]+ (23) 



D 2+ (N,p) = (E 2 +_ p + E° p ) - E 2 N + , (24) 



respectively. However, in the fission processes (z\ > 0) as in Eq.(||), a negative DE does not mean that the cluster 
would decay. It is because of the existence of a fission barrier which originates from the short-range attractive (due 
to binding energy) and long-range repulsive (due to coulomb repulsion) forces between the charged products. The 
situation is shown in Fig. |lj. Any excitation above the barrier, which may be induced by collisions or radiation, will 
eventually make the expected decay possible. One of the main deficiencies of the JM is that it gives negative valuestl 
of surface energies for r s < 2. Our using of SJM-SC is expected, therefore, to give more realistic values of surface 
energies and barrier heights. The fission barrier B% (N,p) is approximated by the Coulomb interaction E c of two 
touching spheres (i.e., the fission products) and the DE as 

B*(N,p) = DZ 1 (N,p)+E c (25) 

For the Coulomb interaction between the two fission products we take into account their polarizabilities. n The 
interaction energy of two charged conducting spheres can be calculated numerically using image charge methodE An 
equally good but much simpler approach is the use of the analytical expressionr 9 IP°l for the interaction between charges 
z\ and Z2 with polarizabilities a\ and a 2 at a distance s 

nn i , , Ai 1z\ zi a\ 1 Zn an 1 z? 

^ = ^-^7-^7- (26) 



Here Aj is given by 



4 = 1-^, (27) 



and the polarizahility of a conducting sphere (i.e., the metal cluster) with radius R i E a = R 3 . An other formula 
which was useduS for the Coulomb interaction of two touching conducting spheres is given by 

El2 = 2ziZ2 - (28) 

' Ri + R 2 + 2SR v ; 

where, for silver, the value 6R = 0.94 takes the polarizability into account. The BH's for small clusters —/obtained 



from this formula, are somewhat smaller than those we obtained using Eqs. (P6|), (27). Koizumi et. aioo have 
calculated the barrier heights for the fission of doubly charged silver clusters using a shape function in the LDM with 
shell correction. 



III. RESULTS AND DISCUSSION 



After an extensive self-consistent SJM-SC calculations, we have calculated the equilibrium r s values and the energies 
of Agjy (2=0,1,2,3,4) for different cluster sizes (1 < N < 100). To show the main differences in the equilibrium r s 
values of these clusters, which are appreciable for relatively small clusters, we have plotted, in Fig. ||, the corresponding 
f 8 (N) values only up to N = 34. As is obviously seen in the figure, the neutral and singly ionized clusters are self- 
compressed for all values of N. This is because of the dominant effect of the surface tension. However, for multiply 
charged clusters, the f s (N) values cross the bulk border (i.e., r s = 3.02) at some N which we show it by AT . Our 
results show that, in general, for larger values of charging, the self-expansion persists up to larger values of No. That 
is, N$ + < < N£ + < ■ ■ For example, here we have obtained the values of 7, 17, 23 for N 2+ , , N* + , 

respectively. This means that, for larger charging values, the coulomb repulsion between excess charges dominates 
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the surface tension up to larger values of Nq. It is also clearly seen in Fig. || that for clusters with the same numbers 
of electrons but different numbers of atoms, N, the following inequality holds 

f° s (N) < fl+(N + 1) < f 2 +(N + 2) < f 3 +(N + 3) < • • • . (29) 

For Ag^ 1 " clusters, we could not find any solution of the Eq. (^lj) for N — 5. That is, the single remaining electron is 
not able to bind the 5 constituent ions to each other in the Agg system. However, the solutions of ^ + (6) = 14.76 



and "rf + (4) = 29.0 have been obtained for Eq. (21) which are so large that one can not realize the corresponding 
bound states experimentally and we ignore these bound states. 

Figure || shows the equilibrium energies per atom in electron-volts for Ag%, Agjy", Ag^" , Ag^, and Ag^ with 
different cluster sizes (1 < N < 34). For comparison, we have also plotted the bulk value (e = — 7.89eT^) by a dashed 
line. As is seen, by increasing the charge of a given iV-atom cluster, the coulomb repulsion between the excess charge 
induces an inflation in the cluster (see Fig. ^) and therefore, the density of the material in the cluster decreases which, 
in turn, leads to a smaller binding energies. 

In Fig. [|(a) we have plotted the DE's of the most favored decay channels for the process Ag]~}~ — > Ag]^_ + Ag^. 
By definition, the DE is minimum in the most favored channel. We have shown the most favored value of p by p* . 
The solid small square symbols show the most favored values p* on the right vertical axis whereas, the corresponding 
DE's, Dq + (N,p*), are shown on the left vertical axis by large open squares. The dashed line is the result of a fitting 
to the quantal DE's. As is seen, the magority of the clusters have positive DE's and therefore, they are stable against 
the spontaneous decay. However, the remaining clusters have negative DE's and accorgingly, they decay into smaller 
fragments. Clusters close to the closed-shell ones, decay by emitting a monomer or dimer. On the other hand, clusters 
that are far from being a closed-shell, can break into two fragments each of which are close or identical to closed-shell 
ones. For example, AgJ + emits a neutral monomer and the remaining is a singly charged dimer, Ag^ emits an Agg 
which is a closed-shell and the remaining is Ag^ which is close to a closed-shell, and finally, AggJ emits Ag°Q and 
the situation is similar to the latter one. Except for the closed-shell singly ionized cluster Agg^J", all other closed-shell 
singly ionized clusters, Ag^ (N—3, 9, 19, 21, 35, 41, 59, 91, 93) are stable against the spontaneous decay. The dashed 
fitted line which resembles the result of liquid-drop model (LDM) calculations ( see Fig. 4 of Ref. jL{|), predicts that 
all singly ionized clusters are stable and trpnasymptotic value of DE is constant and equal to 0.20eV. 

Figure [|(b) compares the experimental^ dissociation energies with the monomer DE's D^ + (N, 1), dimer DE's 
Dl + (N, 2), and the most favored DE's Dl + {N 7 p*). As is seen, the most favored fragments are somewhere monomers, 
somewhere dimers and somewhere none of them. The general trend of the calculated monomer dissociation energies 
is similar to the experimental one and has a better agreement with the experiment than the other two DE's. That 
is, from N — 3 to N = 4 the energy decreases; from N = 4 to N = 9 the energy increases in the mean; a decrease on 
going from 9 to 10; an increase from 10 to 21; and finally, a decrease from N — 21 to N — 22 and again increasing 
from N — 22. However, our results lack the odd-even staggering because, it originates from_the nonspherical shapes 
for the jellium. Resorting to non-spherical shapes also decreases the pronounced shell effects.Ej The relative smallness 
of our calculated DE's can be explained in terms of the details of the experimental setup. If the experiment starts 
using neutral Agjy clusters, then the equilibrium r s values would be smaller than the bulk value (see the Z — plot 
in Fig. |2|), and therefore, the total energies would be more negative (see the Z = plot in Fig. ||). Now, irradiating 
the parent neutral cluster with a high power laser beam would lead to the ionization of the neutral cluster. If the 
photons also interact with the ionized cluster before the ionized cluster achieves its relaxed state, then the equilibrium 
r s value of the ionized cluster would be less than the relaxed value (In our calculations we have used the relaxed 
values at all steps.), and therefore, the magnitude of the total energy of the ionized parent cluster would be larger. 
This fact would lead to larger values of the DE's. Comparing the experimental data with our results we conclude that 
the photo-dissociation occurs before the relaxation of the parent ionized cluster is completed. The other extreme is 
that we consider the 1 sudden'' approximation in which we assume that the relaxation time for the ionized cluster is 
infinite and the cluster undergoes the dissociation without changing the volume (i.e., the saturation approximation 
which is used in nuclear fragmentations and ordinary jellium calculations for clusters). In reality, neither of these 
extreme 1 relaxed or 'sudden? approximations are at work but something in between. 

In Fig. ||(a) we have shown the most favored products Ag°* and the dissociation energies Dq + (N,p*) for the 
decay of Ag^ via evaporation channel. For this process, as in singly ionized case, the dashed fitted line predicts no 
spontaneous decays and shows a higher constant asymptotic DE as 0.40ey. The most favored products are mainly 
monomers, dimers and octamers. 

Figure ||(b) shows the barrier heights Bf+(N,p*) for the most favored channels of the process Ag^ — » Ag]^_ p + 
Agp + . By definition, the most favored fission channel has a minimum value for the BH. As is seen, some of the BH's 
are negative. The negativity of a BH means that we need no energy to supply the system to initiate the fission. The 
dashed line which shows the mean behavior of BH intersects the zero line at j\P nean pa 29 (the mean appearance size) . 
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This means that on the average, all Ag^ clusters with N < 29 are unstable against spontaneous fission. The values 
of p* show that most of the emitted fragments are closed-shell Ag]^~ clusters with N = 3, 9, 21. 

In Fig. 0(c) we have compared the most favored decays of Figs. 0(a) and 0(b). It is clearly seen that in a certain 
size range, the fission and evaporation definitly start their competition. Our quantal results in Fig. 0(c) show that in 
the size range 21 < N fjr-26 the evaporation dominates the fission which is in good agreement with the Katakuse et. 
al. experimental resultsEj that reveal fission for N < 22. However, our result is slightly larger than the Kruckeberg 
et. alt^i experimental data which show that the fission occurs for N < 16. This difference in the experimental results 
depends on the details of the experiment. For N > 26, our results show that fission dominates again. To estimate 
the size range at which evaporation completely dominates the fission, we simply find the intersection point of the 
two mean behaviors (dashed lines) in Figs. 0(a) and 0(b). A simple calculation gives this mean critical value as 
jV™ » 50. That is, in an induced fragmentation experiment of Ag^*~ clusters, the evaporation dominates the fission 
for N > 50. 

In FigJ|(d) we have compared the most favored values Dq + (N, p* ) and B 1 1 (N, p* ) with the experimental threshold 
energiesLa Here, also we have smaller DE's and BH's compared to the experiment. One reason for this behavior is 
that the equilibrium volume of the parent cluster is not equal to the sum of the equilibrium volumes of the product 
clusters (i.e., the 1 relaxed approximation) but is larger. The larger value of the equilibrium r s leads to a smaller 
magnitude of the initial energy and therefore, by Eq. (|2^) to a smaller barrier heights. In other words, the energy 
needed to deform the parent cluster toward the fission ( In deformation the surface area increases.) is partly paied as 
a result of self-expansion of the parent cluster. Our calculations show that in almost all decay channels, the sum of 
volumes of the decay products is smaller than that of the parent cluster which can be explained by the fact that in 
smaller clusters the surface effect is higher than that in larger clusters. The calculated results for the most favored 
values of DE's at N = 9 and N = 10 are very close to the experimental values. The most favored products at N = 9 
and N = 10 are neutral dimer and monomer, respectively (see Fig. 0(a)). 

In Fig. 0(e) we have compared the calculated monomer DE's Dq + (N,1) and the singly charged trimer BH's 
Bl+(N,3) with the experiment. Here also the difference is appreciable. 

In Fig. 0(a) we have plotted the p* and the Dq + (N,p*) for different cluster sizes. The situation is similar to other 
previous evaporation processes. Here also the asymptotic behavior of the fitted line predicts no decay and has a 
constant value of about 0.50eV\ 

In Fig. 0(b), we have plotted the BH's Bf+(N,p*) for the most favored channel of the binary fission of the process 
Ag 3 j^~ — > AgJ^L p + Agp + . Here, the fission products with smaller charge are more or less the same as those in the 
fission of Ag^*" clusters. The mean behavior dashed line intersects the zero axis at _/V™ can rj 53. The slope of this line 
is larger than that of Fig. 0(b). 

Figure 0(c) compares the most favored decays of Figs. 0(a) and 0(b). To our knowledge, there is no experimental 
results in the literature on the decay of Ag^ with Z > 3. It is seen that at N — 38, evaporation dominates and from 
N = 39 to N — 42 evaporation and fission are equally probable. From N = 43 to N = 58, except for N = 54, fission 
dominates again. From N = 59 to N = 67 the evaporation process overcomes and so on. To obtain the mean critical 
value, we find the intersection point of the two mean behaviors (dashed lines) in Figs, 0(a) and 0(b) which results in 
the value N™ an w 80. 

Fig. ^(a) plots the p* and the D + (N,p*) for the evaporation processes of Ag^. In this figure, one notes the 
unstability of two smallest sized Agy + , Agg + and the stability of almost all others against the evaporation. The fitted 
dashed line predicts no evaporation and has the asymptotic value of 0.75ey. The evaporation products are seen to 
be mostly neutral dimers and a few monomers and octamers. 

In Fig. ^(b) we have plotted the p* and the Bf+(N,p*) as functions of N. We see that at N = 21 the BH becomes 
positive for the first time. The dashed line shows the mean behavior of the fission barries. This fitted line has crossed 
the zero axis at 7V{ ncan rs 74. That is, on the average, the Ag^ 1 " clusters are stable against the charge-asymmetric 
fission channel for N > 74. The most favored charged products are mostly magic clusters, Ag]^~ with iV = 3, 9, 21. 

The results of charge-symmetric binary fission of Ag^ for most favored decays are shown in Fig. 0(c). The smallest 
positive BH occurs at N — 22. The mean behavior (the dashed line) intersects the zero line at N™ oan w 51. 

In Fig. 0(d) we have compared the results for the three different decay processes of Ag^. It is seen that these 
clusters smaller than N ~ 50 are highly unstable. As is seen from the figure, it is difficult to specify the competitions 
for the quantal values. However, to give quantitative values, we have compared the fitted lines in Fig. 0(e). It is seen 
that for N < 20 the charge-symmetric fission is the dominant spontaneous decay process but, for 20 < N < 51 the 
dominant spontaneous fission process changes to the charge-asymmetric one. For 51 < N < 73 the only spontaneous 
fission decay process is the charge-asymmetric one. Clusters larger than jV™ ean ks 73 are stable against any spontaneous 
decay process. In an induced fragmentation experiment of Ag^ 1 " clusters, the dominant process for 73 < N < 107 is 
charge-asymmetric fission, and for N larger than jV™ ean = 107 the evaporation process dominates. To summarize, for 
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smaller clusters the charge-symmetric fission is dominant, and larger clusters prefer to decay via a charge-asymmetric 
fission process. 

Besides the most favored quantities which are strongly related to the stability of the charged cluster and were 
explained in the above lines, it is also interesting to calculate the DE's and BH's for a process in which the fragment 
products are specified. Consider the process Agfj — ► Agf ! _ 1 + Ag° in which one of the products is a neutral monomer. 
We have calculated the dissociation energies Dq(N,1) for all values of Z=l, 2, 3, 4 and N < 100. The calculated 
values show pronounced shell effects as in previous figures of the most favored channels. However, the mean behaviors 
have asymptotic constant values. For Z = 1, 2, 3, 4, these asymptotic values in electron-volts are 0.95, 1.00, 1.08, 1.15, 
respectively. This means that the monomer evaporation from a singly charged cluster needs a smaller energy than 
from a doubly charged and so on. The same analysis for the dimer evaporation in the process Agj^ — -> Ag^_ 2 + Ag° 
shows also constant asymptotic mean behaviors for the Dq(N, 2). The obtained values in electron-volts are 0.49, 0.59, 
0.72, 0.87 for Z—l, 2, 3, 4, respectively. This means that, as in the monomer evaporation, the detachment of a dimer 
from singly ionized cluster is easier than from a doubly ionized cluster and so on. However, comparing the dissociation 
energies for monomer and dimer evaporation (keeping the charge constant) shows that atomic evaporation needs more 
energy than dimer evaporation. 

Now, we consider the fission processes Agfj — > AgfjZ^ + Ag 2 + and Ag^- — > Agj^T 3 + AgJ + in which one of the fission 
products is a singly ionized dimer or a singly ionized trimer. The mean behaviors of the BH's for these processes are 
plotted in Fig. [3| As is seen, the energy needed to detach a singly ionized dimer decreases by increasing the charge 
of the parent cluster. This behavior should be contrasted to the behavior in the monomer or dimer evaporations. We 
recall that in the monomer or dimer evaporation, the dissociation energy increases by increasing the charge of the 
parent cluster. 

It is now easy to find the mean sizes at which atomic evaporation process dominates the fission into singly ion- 
ized dimer or trimer for each charging value of the parent cluster. In doubly charged silver clusters, the monomer 
evaporation dominates the singly charged dimer and trimer detachments at N = 11 and N — 31, respectively. The 
corresponding numbers for triply charged clusters are 21 and 66. For parent clusters Ag^*~, the numbers N = 38 and 
N = 120 have been obtained. To summarize, by increasing the charge of the parent cluster the competition occurs at 
larger values of N. 

IV. CONCLUSION 

In this work, we have studied the fragmentations of multiply charged silver clusters taking into account the structural 
relaxations of the neutral and charged parent as well as daughter clusters. To calculate the relaxed-state sizes and 
energies of the clusters we have employed the stabilized jellium model with self-compression using a spherical geometry 
for the jellium background. Using these relaxed-state radius and energy for the clusters, we have calculated the 
dissociation energies and barrier heights for evaporation and fission processes in all possible channels. For the barrier 
heights, we have used the two-touching-spheres model with taking into account the polarizabilities of the two charged 
products. Comparison of our most favored results with the experimental data shows that our results lie under the 
experimetal results but, the critical size for the competition of the evaporation and the fission of doubly charged silver 
clusters is predicted in good agreement with the experiment. This comparison also reveals that the fragmentation 
processes mostly occur before the complete relaxation of the charged parent clusters. That is, in the above-mentioned 
experiments the structural relaxation time is larger than the average time elapsed for the fragmentation of the ionized 
parent cluster. Having the initial (just after ionization) and the relaxed sizes rf ini , rf lela of a Z-ply ionized cluster, 
one may choose an rf frag value (rf ini < rf frag < rf rela ) for the ionized cluster (just before the fragmentation) such 
that the calculated values coincides the experimental ones. Then using a linear interpolation it is possible to calculate 
the relative fragmentation time for an ionized cluster. In ordinary jellium model calculations, the assumption 

r° s (N) = rf(N) = rf ini (N) = rf^ g (N) = rf Icla (N) = r s , bulk (30) 

is used. It should be mentioned that for exact matching of the calculated and experimental values one should use 
non-spherical shapes. 

We have obtained the asymptotic DE's for the most favored channels in evaporation processes by fitting a simple 
curve on the quantal results. The result shows that the asymptotic values increase by increasing the charge of the 
parent cluster. In the case of Ag^*~, we have shown that for relatively small clusters the charge-symmetric fission 
process is dominant and then, before dominating the evaporation process the charge-asymmetric fission process 
overcomes. In general, the critical size (at which the evaporation dominates the fission) increases by increasing the 
charge of the parent cluster. The results show that the neutral p-mer dissociation energy increases by increasing 
the charge of the parent cluster; and for a given charged parent, the atomic evaporation needs more energy than a 
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dimer evaporation. Finally, it has been shown that the energy needed for the detachment of a singly charged dimer 
or singly charged trimer decreases by increasing the charge of the parent cluster. However, for a given parent cluster, 
the detachment of a singly charged trimer is easier than that of a singly charged dimer. 
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FIG. 1. Fission barrier in the two-touching-spheres model. B, D, and E c are barrier height, dissociation energy, and the 
coulomb energy between two touching charged conductors, respectively. 

FIG. 2. The equilibrium r a values in atomic units for Z-charged clusters as a function of the cluster size N. The dashed line 
is the bulk Wigner-Seitz radius for Ag (r s = 3.02). 



FIG. 3. The equilibrium total energies per atom in electron volts for Z-charged clusters as a function of the cluster size N. 
The dashed line is the bulk SJM value for Ag (e = -7.89eV). 

FIG. 4. (a)- The right vertical axis shows the size of the fragment in the most favored channel, p* , of singly ionized clusters 
by small solid squares. The left vertical axis shows the dissociation energies, DE, in electron volts for the most favored decay 
channels of singly ionized clusters by large open square symbols. The dashed line is the fitting result which has the asymptotic 
value of 0.20eV. (b)- The dissociation energies for monomer (p = 1), dimer (p — 2), and the most favored (p = p*) p*-mers 
are compared with the experiment. The figures on the experimental data points specify the fragment size p. The experiment 
shows an odd-even staggering. 

FIG. 5. (a)- The same as in Fig. [l|(a) for doubly ionized clusters, with the asymptotic value of the dashed fitted line as 
OAOeV. (b)- The most favored size of the fragments, p* , and the barrier heights for this most favored values in electron volts as 
function of the cluster size N. The dashed fitted line intersects the zero axis at N — 29. (c)- Comparison of the most favored 
values of the DE's and the BH's. is Dg + for Z\ = and B^t for z\ = 1. Evaporation starts from N = 21. (d)- Comparison 
of the most favored DE's and BH's with the experiment. The figures on the experimental data points show the fragment sizes, 
(e)- Comparison of the monomer DE's and the singly ionized trimer BH's with the experiment. 

FIG. 6. (a)- The same as in Fig. ^(a) for triply ionized clusters, with the asymptotic value of the dashed fitted line as 
0.50eV\ (b)- The same as in Fig. |fj(b) with the intersection point at N = 53. (c)- The same as in Fig. ^(c) with starting the 
evaporation from N = 38. 



FIG. 7. (a)- The same as in Fig. ^(a) for 4-ply ionized clusters, with the asymptotic value of the dashed fitted line as 0.75eV. 
fb)- The same as in Fig. |^(b) for charge-asymmetric fission with the intersection point at Ni = 74. (c)- The same as in Fig. 
pKb) for charge-symmetric fission with the intersection point at N2 = 51. (d)- Comparison of the evaporation DE's with the 
charge-asymmetric (zi = 1) and the charge-symmetric (zi = 2) fission BH's. (e)- Comparison of the fitted lines for the processes 
in (d). At N = 20, the charge-symmetric and the charge-asymmetric fissions compete. For N > 51, charge-symmetric and for 
N > 73, the charge-asymmetric fissions stop, respectively. For JV™ ean > 107 evaporation overcomes the fission processes. 

FIG. 8. Comparison of the fitted values of the BH's for dissociating a singly ionized dimer or trimer from a Z-p\y charged 
cluster. 
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